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Abstract 

We calculate the dynamical relativistic corrections, originating from radiative one- 
gluon-exchange, to the leptonic decay width of heavy quarkonia in the framework of a 
covariant formulation of Light-Front Dynamics. Comparison with the non-relativistic 
calculations of the leptonic decay width of J= 1 charmonium and bottomonium S-ground 
states shows that relativistic corrections are large. Most importantly, the calculation of 
these dynamical relativistic corrections legitimate a perturbative expansion in ag, even in 
the charmonium sector. This is in contrast with the ongoing belief based on calculations 
in the non-relativistic limit. Consequences for the ability of several phenomenological 
potential to describe these decays are drawn. 
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1 Introduction 


Although the structure of heavy quarkonia in terms of a heavy quark-antiquark non-relativistic 
bound state is known for a long time, many recent developments show the importance of a 
relativistic description of such states. Mainly two paths are followed: 

i) One can first estimate relativistic corrections to the qq bound state in a n/c expansion. 
This is the aim of the NRQCD formalism |^. When the quark mass is large enough, this 
is certainly adequate since this expansion should converge rapidly. It is however not clear 
whether the charm quark mass is heavy enough for such an expansion to be valid. To 
check it out, one should compare the results of NRQCD with relativistic calculations. 

ii) Relativistic corrections can also be calculated in a relativistic framework to describe the 
two-body bound state, such as the Bethe-Salpeter formalism, or Light-Front Dynamics 
(LFD). In a previous article |0], we have investigated the relativistic kinematical correc¬ 
tions to the leptonic decay width of the J/'ip and the T induced by the finite momentum 
between their quark and antiquark. We used an appropriate formalism, the so called 
’’Covariant formulation of Light-Front Dynamics” (CLFD), the details of which can be 
found in Ref. [§. As pointed out in this reference, the CLFD allows a straightforward 
comparison between relativistic and non-relativistic calculations. 

We shall investigate in this article the importance of dynamical relativistic corrections. Kine¬ 
matical relativistic corrections can be calculated without any knowledge of the dynamical ori¬ 
gin of the two-body wave function. Dynamical corrections, on the other hand, correspond to 
relativistic corrections to the wave function itself. While the non-relativistic two-body wave 
function of the J = 1 state has only two components (the S and D states), the relativistic wave 
function formulated on the light-front has 6 dynamical components 0 . If the dynamical origin 
of the two-body interaction is known, these components can, in principle, be calculated exactly 
or else calculated in perturbation theory, starting from the non-relativistic components. This 
has been done for instance in Ref. IQ, for the calculation of the relativistic corrections to the 
deuteron wave function. 

In the case of heavy quarkonium, the two body interaction has two main parts: a conhning 
interaction, and the interaction coming from the exchange of an effective gluon (One-Gluon- 
Exchange, OGE). The latter interaction gives rise, in the non-relativistic limit, to the well known 
Goulomb interaction. While little is known on the dynamical origin of the conhning potential, 
one may hope that for heavy quarkonia, the two-body wave function is mostly sensitive, for 
relativistic corrections, to the short range part of the interaction, i.e. to the one-gluon exchange 
potential and its 1/r non-relativistic behaviour. 

This is the usual assumption made when calculating the dynamical relativistic correction 
coming from the OGE interaction. In zeroth order in a k‘^/m? expansion, where k is the 
relative momentum of the quark-antiquark pair in the two-body bound state, it gives rise to 
the following well known correction to the leptonic decay width of J = 1 states P, |^: 

ry = rr (i - ^) . (1) 
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In the case of the J/'h leptonic decay width, this contribution amounts to a 50% reduction. 
This large correction clearly calls for a relativistic calculation of this contribution to all orders 
in This is the aim of the present article. 

We shall follow very closely the procedure of Ref. p|. We refer the reader to this reference, 
and to Ref. [Q for a review article on CLFD. We introduce in section 2 the scheme of the 
calculation of the leptonic decay width. The dynamical relativistic corrections are introduced 
in section 3. Numerical results are presented in section 4, and our conclusions and perspectives 
are drawn out in section 5. 


2 The leptonic decay width of heavy quarkonium 

2.1 The quarkonium wave function in CLFD 

Among the various ways to deal with relativity in the description of bound states, we will 
focus on LTD. In the standard formulation of LTD 0, the wave function of a bound system 
is defined on a plane characterised by the equation t + z/c = 0. The usual Schrodinger, equal 
time, formalism is easily recovered in LTD by letting c go to inhnity. 

The description of relativistic heavy quark systems in LTD has many advantages, of which 
the most important one may be the absence of vacuum fluctuations. Accordingly, a meaningful 
decomposition of the state vector describing the system in terms of Fock components of dehnite 
number of particles is possible. Of course, the number of Fock components to be considered 
in any practical calculation depends on the dynamics of the system, and on the kinematical 
regime one is interested in. 

The most serious drawback of this formulation, however, is that the position of the light- 
front t -|- z = 0 (with c = 1) is not rotational invariant. Since rotations in the zx and zy planes 
change the position of the light-front, the associated generators shall depend on the dynamics 
and cannot be reduced to kinematical transformations [§]. This means that one needs to know 
the complete dynamics in order to write down the general structure of a bound state of definite 
angular momentum. It also means that any electromagnetic operator should have the same 
(dynamical) transformation properties to match the bound state wave function one’s. This is 
essential to guarantee that any physical amplitude (or cross-section) does not depend on the 
particular choice of the light-front we start with. 

We therefore need an explicit procedure to exhibit in a convenient way these dynamical 
transformations. This is achieved in the covariant formulation of LFD. Our starting point is 
the invariant definition of the light-front by cn-x = 0 where a; is a (unspecified) light-like four 
vector (o;^ = 0). The standard formulation of LFD can be easily recovered with the particular 
choice uj = (1, 0, 0, —1) for the light-like four vector. 

This definition of the light-front is explicitly invariant by any four-dimensional rotations, 
or any three-dimensional rotations and Lorentz boosts. Consequently, these transformations 
become cu-dependent, but do not necessitate the knowledge of the dynamics of the system 
to construct them explicitly. The dynamics enters now into the cu-dependence of the wave 
function and of the electromagnetic operator, in such a way that any physical amplitude should 
not depend on the particular position of the light-front, i.e. on u, unless approximations have 
been made. In this case, which is almost always true in practice, the explicit covariance of the 
approach enables us to exhibit the cn-dependence of the amplitude and to separate the physical 
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part from the non-physical, a;-dependent one, as we shall explain below for the leptonic decay 
amplitnde. 

The wave fnnction $, of a two-body bonnd state, can be decomposed in terms of all 
the possible independent spin strnctnres compatible with the qnantnm numbers of the state 
studied. In the particular case of vector mesons, which we are interested in, we can write : 

= y/mel{p)%f^{k2Wv^^{ki) , (2) 


with: 
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This decomposition is similar to the decomposition of the deuteron wave function found in 
Ref. IQ. The six components of the wave function, depend on two invariants. In order 

to make a close connection to the non-relativistic case, it will be convenient to introduce another 
pair of variables 0 dehned by: 
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L-\V)ki = h- 
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(4) 
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(5) 


where: 


V = ki + k2 . (6) 

The relativistic momentum k corresponds, in the frame where ki + k 2 = 0, to the usual relative 
momentum between the two particles. Note that this choice of variable does not assume that 
we restrict ourselves to this particular frame. The unit vector n corresponds, in this frame, to 
the spatial direction of uj. In terms of these variables, the wave function takes a form similar 
to the non-relativistic one: 




(7) 
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where Wo- is the two-component Pauli spinor normalised to = 1, and a are the usual Pauli 

matrices. The relation between and is the same as the relation between the spherical 
function Y^{n) and n. 

The coefficients of the spin structures in Eq. @ and Eq. (H) are scalar functions of two 
independent invariants, which we can choose to be k'^ and k-n, since these variables are only 
rotated by a Lorentz boost 0]. In the non-relativistic limit, only two components remain, /i 
and / 2 , and they only depend on k"^. This can be easily seen if we keep track of the c factors, and 
then let c goes to inhnity to get the non-relativistic limit. In this study of heavy quarkonium 
states, one may neglect the tensor component, / 2 , so one is left with the non-relativistic wave 
function /i = The relation between (pi, (p 2 and is given by |^: 
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where Ek = VP + m?. Note that the wave function /i = is normalised according to: 

2 dPk m 






= 1 . 


( 11 ) 


2.2 The leptonic decay width 

Since the total physical amplitude for the process Mj e+e“ can be factorised into two 
separate amplitudes Mj —> 7 and 7 —*• the relevant physical information is completely 

included in the amplitude to produce a photon with polarisation from a vector quarko¬ 
nium state of polarisation e^. 

Since our formulation of LED is explicitly covariant, we can decompose in terms of 
all possible tensor structures build up with the two four-momenta at our disposal, p and cu, 
where p is the four-momentum of the quarkonium P|. Thus, we can write: 


with: 


A^w = Fapp + 


Bi 

2uj-p 




B 2 
2a;-p 


af + B, 




{uj-p)' 




( 12 ) 


^lp 

'1 

pp P^F 

y M2 ’ 

(13a) 

pp 

'2 

= p^Uj’’ + P^UJ^ , 

(13b) 

pp 

'3 

= p>^ujp - pPujf^ , 

(13c) 

PP 

'4 


(13d) 

pp 

'5 

p^pp 

M2 

(13e) 


We have denoted by M the mass of the quarkonium. Note that the term proportional to D 
does not contribute to the leptonic decay width. The amplitude has two kinds of terms. 
The hrst one, proportional to F, is the physical contribution to the decay width. The other 
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Figure 1: Leading order contribution to the leptonic decay width. The dashed line represents 
the spurion (see text for details). 


three, proportional to i?i, B 2 and are a;-dependent contributions. In an exact calculation, 
the coefficients i?i, B 2 and B^ should be zero since the physical leptonic decay width should 
not depend on the particular orientation of the light-front one starts with. However, in any 
approximate calculation, these terms may be non-zero, but they are not physical. Consequently, 
they have to be eliminated in the computation of the physical leptonic width. 

To extract F from the general amplitude one can hrst multiply successively 

by the five tensor structures Oi to 05 given in Eqs. ( 0 . This gives a system of five coupled 
equations which is solved to get the physical amplitude F. Thus, we hnd: 

F = 1 (Ji - 2 I 2 + h + h) , (14) 

with: 


h 


At 


pp 


(15) 


These quantities are easily evaluated once Ai^^ is calculated from the process under consid¬ 
eration, using the diagrammatical rules given in Ref. |^. Given F, the decay width can be 
calculated [ 0 , and is given by: 


T 


dvr 

3ii^ 



(16) 


where Cg is the electric charge of the quark and a is the electromagnetic fine structure constant. 


2.3 The zeroth order calculation 

The leading order contribution to the leptonic decay is shown in Fig. where we have removed 
for simplicity the trivial vertex 7 —>■ Off-energy shell effects are governed by the variable 

r, which is unambiguously determined by four-momentum conservation and the on mass-shell 
condition for each particle [^. Using the conservation law at the bound state vertex, we haveQ 

ki + k2=p + 0JT . (17) 

To keep track of this conservation law, we represent the four-vector ujt by a dashed line in 
every diagram (the so-called spurion line, see Ref. [Q for more details). Note that the outgoing 

^One should note that in CLFD, in the frame where we have ki -\- k 2 = 0, the total momentum p is not 0. 
Instead p + wr is 0. 
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photon has virtuality = p^. However, it can be assumed to be a ’’physical particle” in 
the sense that the hnal process 7 —*• e+e“ is completely disconnected from the decay process 
Mj —> 7 , and is exactly calculable in terms of outgoing free electron and positron. As already 
explained in details in Ref. [@], the zeroth order contribution to the amplitude, represented in 
Fig. 0, is given by: 

/ (]3u 

———----Tr [7^(^i+ m)0^(m - 1 ^ 2 )] , (18) 

(27r)'^2£fci(l - X) 


where x = oj-ki/oj-p and (f)'^ is defined in Eq. (|^). The physical part of the amplitude in zeroth 
order calculated by Eq. ( 0 , Fq, can now be written in the form: 


^^0 = 


(Fk m 
{27i)^ek 


Oo{P)4>^\k^) , 


(19) 


where Oq(JF) is: 


Oo{P) 




2 m 




( 20 ) 


Note that in leading order in a k‘^/nF expansion, kinematical relativistic corrections originate 
only from the factor m/sk of the relativistic phase space volume. In the non-relativistic limit, 
one has: 


= —2\/6m(j)^^{r = 0 ) . 


( 21 ) 


3 Dynamical relativistic corrections 

3.1 Radiative corrections 


As already pointed out, the dynamical relativistic corrections necessitate the knowledge of 
the dynamical origin of the two-body wave function, he., the way quark-antiquark states are 
bound. This goal is still far from being achieved. The standard assumption to overcome 
this issue is to suppose that for quark masses heavy enough, the dynamics is governed by a 
perturbative one-gluon-exchange. The physical amplitude for the relevant processes are 
indicated in Fig. ^ They can be calculated analogously to Eq. ([T8|) using the diagrammatical 
rules given in Ref. [§]. The details of the calculation are too lengthy to be shown here, but 
present no particular difficulties. The amplitude is both ultraviolet and infrared divergent. The 
computation of these contributions is detailed in appendix. 

Since these diagrams are ultraviolet divergent, we also need to include the renormalisation 
of the quark charge at the electromagnetic vertex, as has been already detailed in Ref. 0 . The 
corresponding diagram for this contribution is shown in Fig. |^. The counter-term Z is given 
by: 


Z 



9 1 / 

^ + 73- log — 

STT 271 \m^ 





( 22 ) 
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Figure 2: Radiative correction to the leptonic decay width. 



Figure 3: Vertex renormalisation contribution to the leptonic decay width. 


where the leading | is a factor of colour generated by the gluon exchange. 

The ultraviolet regularisation of the amplitude is done according to the Pauli-Villars 
prescription, with a cut-off A. This is the same prescription as in Ref. |^. The amplitude is 
made infrared finite by giving a small non zero mass fi to the photon. We shall see in the next 
section that the hnal contribution remains infrared hnite, as it should. 

3.2 Coulomb subtraction 

As it is well known, one should be very careful in the evaluation of dynamical relativistic 
corrections as indicated in Fig. H. Indeed, when we evaluate the zeroth order contribution 
shown in Fig. |^, with the non-relativistic wave function, solution of the Schrodinger 

equation with the non-relativistic two-body potential: 

= ^oge{^) + Vconfir) , (23) 

we implicitly incorporate corrections of the type shown in Fig. ^ in leading non-relativistic 
order at least. One should therefore make sure that these contributions are properly removed 
from the relativistic calculation in order to avoid double counting. 

The usual procedure to do this is to remove from the relativistic OGE the Coulomb 
interaction. We shall see below that our formulation of the leptonic decay width in CLFD 
enables us to have a clear handle on the contribution to remove, at any level of approximation. 
Let us start from the Schrodinger equation, written in the form: 

- -7^ j . (24) 

with Vqq given by Eq. (^3|), and s = 4(A;^-|-m^). The contribution to which corresponds 

in perturbation theory (hrst order in a*) to the OGE interaction in the non-relativistic limit is 
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Figure 4: Coulomb contribution to the radiative correction of the leptonic decay width. 


therefore given by: 




4 4m 


d^k' 


a.. 




3 s - M2 7 (27r)3 _ gy 


Inserted in the zeroth order calculation indicated in Fig. |I|, and written in Eq. 
contribution, 6F, to the physical amplitude: 


(25) 

, it gives the 


SF^ f , (26) 

J [ZTiyek 

This contribution is represented graphically in Fig. It is ultraviolet hnite, but diverges in 
the infrared region. We regularise the latter divergence by giving a small hnite mass /i to the 
gluon. Note also that in Eq. (^), the phase space volume in k' is the non-relativistic one since 
it involves the non-relativistic wave function, solution of the Schrodinger equation. However, 
in Eq. (^), it involves the relativistic phase space volume since one has to perform the integral 
over k on the whole momentum range, as done in Eq. (^) for the calculation of kinematical 
relativistic corrections. 

The contribution dF of Eq. (^61) , together with 6(j)^^{k‘^) given in Eq. (^5[) , should be 

removed from the total relativistic amplitude given by Eq. (O). In the non-relativistic 

-^2 

limit for the two-body bound state (limit [k' /nF) -C 1), it gives the following contribution to 
F: 

o ^ 


3.3 Non-relativistic limit 


In the non-relativistic limit for the two-body bonnd state, the contribution of Fig. to the 
physical amplitude is: 


rNR _ 
^OGE — 


m 

LP 


7 1 , 

- \ -log 

Stt 27r ^ 


( 4 ) 

\m^ J 


1 , (s? 


m^ 


:^NR 


(28) 


while the contribution of Fig. 0 gives, from Eq. 
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Adding the contributions of Fqq^ and then removing the contribution 6F^^, we get the 
hnal, infrared hnite, contribution in the non-relativistic limit: 

p.NR ^ . (30) 

071 

Added to F^^, it gives the well known radiative correction in leading order in as, as indicated 
in Eq. for the total width. 


4 Numerical results 


To easily compare our numerical calculations with the non-relativistic limit found above, we 
consider the following ratio: 


R = 


hi 

rr 


hi To 
To ’ 


(31) 


where Ti is the relativistic decay width to hrst order in ag (calculated in sec. |3.1|) , Tq is the 
relativistic decay width in leading order in (calculated in sec. ^ 
relativistic limit in leading order in ag given by the Van Royen-Weisskopf formula: 


while is the non- 



rr=16^eJ^|0^V = O)r . 

(32) 

We also note: 


(33) 



(34) 


The ratio Rk is simply the kinematical relativistic corrections we already calculated in Ref. [Q]. 
The dynamical corrections have been written in the form of Eq. (^) to explicitly exhibit the 
non-relativistic limit of the radiative corrections when the two-body bound state wave function 
is restricted to very small momenta. In this case, Rjo should go to unity. 

For the hadronic vertex, we use two types of test wave functions to assess the importance 
of these relativistic corrections. On the one hand, we take an harmonic oscillator wave function 
corresponding to a conhning potential for large distances between the quark and the anti-quark. 
The wave function reads: 


= N 



3/4 

exp 



(35) 


In the previous expression, k‘^ is the mean relative momentum squared. In this case, the wave 
function at the origin writes: 


= 0) = N 



(36) 
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On the other hand, we take a Conlomb wave fnnction expected to dominate for very heavy 
quark masses. It is given by: 



As in Eq. (^), is the mean relative momentum squared. Thus, the wave function at the 
origin is: 

= 0) = N^ (klf* . (38) 

Note that the choice of these test wave functions is just a way to investigate the sensitivity of the 
decay constant to various shapes of wave function. It is not intended to imply any assumption 
about the exact dynamics binding the heavy quark-antiquark pair. The normalisation factor 
N is calculated according to Eq. O- It goes to one in the non-relativistic limit. 

In Figs. 1^ and we plot our numerical results for the kinematical and dynamical (ra¬ 
diative) relativistic corrections, for the charmonium and the bottomonium ground state re¬ 
spectively, as a function of the square of the radial wave function at the origin, uq, dehned 
as: 


Mo 


= = 0 ) . 


(39) 


This is indeed the relevant variable in this case since the leading, non-relativistic decay width, 
T^'^, as given in Eq. (|3^ , is directly proportional to it. In the following numerical calculation, 
we choose as = 0.3 for the charmonium sector and as = 0.15 for the bottomonium sector. 

The kinematical relativistic correction is shown in Fig. |^(a) for the J/T, with rric = 


M 


j/'p 


is the same as the one already calculated in Ref. 0 for three different phenomenological 
potentials. The results obtained here for our two kinds of test wave functions are very similar 
in shape. As expected, a larger reduction is obtained for the Coulombic wave function as it 
has larger high momentum components^. The dynamical relativistic corrections are indicated 
in Fig. ^(b). They are much larger, with a sharp decrease as a function of Mq. Again, the 
Coulombic wave function has a larger correction. For a typical value of Uq of the order of 
1 GeV^, the correction is of the order of 0.25. This means that the well known non-relativistic 
radiative correction, — in Eq. (|l|) is reduced by a factor 1/4. This is consistent with the 


hrst estimate of this correction found in Ref. . This result has an important and rather nice 
consequence: it legitimates an expansion in as in the charmonium sector, since the radiative 
correction is now of the order of 10 — 15%, as compared to the 50% reduction given by previous 
non-relativistic calculation with ck^ = 0.3. 

The total correction to the decay width, represented by the ratio R, is shown in Fig. §(c). 
Surprisingly enough, this ratio is rather insensitive to the precise value of Mq to the 

precise value of the radial wave function at the origin. Indeed, the two relativistic corrections 
(kinematical and dynamical) have opposite effects on R and tend to compensate each other, 
leading to an overall reduction factor varying from 0.6 (for a Coulombic wave function) to 


^This correction even diverges in a naive k'^/rn? expansion! 
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Figure 5: Relativistic corrections to the J/tH leptonic decay width, for the two test wave func¬ 
tions. The Coulombic one (solid line) and the Harmonic one (dashed line). These contributions 
are plotted as a function of Mq, dehned in Eq. (|3^ and proportional to the square of the wave 
function at the origin. Sub-hgure (a) gives the kinematical correction, ratio Rk, (b) give the 
dynamical correction, ratio Rd, and (c) show the total correction, ratio R. The mass of the 
charm quark is taken to be 
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in GeV^ 


Figure 6: Same as in Figs. 5 but for the Upsilon decay width. The bottom quark mass is taken 
to Mx/2. The values of Uq for the harmonic oscillator wave function have been deliberately 
restricted in order to correspond to physical ranges for the mean squared momentum between 
the quark-antiquark pair. 
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0.7 (for an Harmonic wave fnnction). This indicates that to have an agreement with the new 
experimental decay width of the J/T HIT 


the sqnare of the radial wave fnnction at the origin 


shonld be of the order of 0.8 — 1 GeV^ 
Refs. 


g, P, [15 


This is the case for several potentials, see for example 
A pnre Conlombic potential at short distance, like the Cornell potential 
is however rnled ont by onr analysis, while a Conlombic potential modified according to 
asymptotic freedom is perfectly viable |g, jl^ . 

The resnlts for the T state are indicated on Figs. j^. The behavionr is very similar to the 
one observed in the charmoninm sector. The overall correction factor is of the order of 0.65 (for a 
Conlombic wave fnnction) and 0.8 (for an Harmonic wave fnnction). The total correction factor 
is snrprisingly large given the large bottom qnark mass. However, the correction factor is mnch 
larger for a Conlombic wave fnnction than an Harmonic oscillator one. This is indeed natnral 
since for large qnark mass the two-body wave fnnction is concentrated to small distances. In this 
region, the high momentnm components of the wave fnnction play a major role in enhancing 
relativistic corrections. 

In order to get an agreement with the experimental decay width, the sqnare of the wave 
fnnction at the origin needs to be of the order of 6 to 8 GeV^. This has several interesting 
conseqnences. On the one hand, two potentials which are compatible with the J/T decay 
width (logarithmic fl^, and power law |T^) give a far too small decay width for the T [j^ 


These potentials have no Conlombic components, and they lack in high momentnm components 
which show np in the T sector. On the other hand, a pnre Conlombic potential still gives a 
too large wave fnnction at the origin. The two potentials of Refs. [^, j^, which correct the 
Conlomb interaction at short distances according to asymptotic freedom, are compatible with 
the experimental T decay widthQ, as well as the J/T one. 


5 Perspectives 

We have presented a coherent formnlation of the relativistic corrections to the leptonic decay 
width of heavy qnarkonia. These corrections inclnde both kinematical and dynamical contri- 
bntions, the latter were calcnlated in first order in ag, bnt to all orders in /w?. 

As already fonnd in onr previons stndy [@], relativistic corrections are very large in the 
charmoninm sector, and not negligible in the bottomoninm sector. Kinematical relativistic 
corrections, which are independent of the dynamics binding the qnark-antiqnark pair, lead 
to a large rednction of the leptonic decay width. Dynamical relativistic corrections, which 
correspond to relativistic corrections to the two-body bonnd state itself, lead to a sizeable 
rednction of the standard correction (—^fy^) fonnd in the non-relativistic limit. This result is 
particularly important since it shows that an expansion in becomes now meaningful. Indeed, 
the first order correction in now amounts to about 10 — 15% correction in the case of the 
J/T, down from 50% in the non relativistic limit. 

While the individual relativistic corrections (kinematical and dynamical) depend strongly 
on the mean momentum squared (or equivalently on the radial wave function at the origin), 
the total relativistic correction is remarkably constant as a function of the square of the wave 
function at the origin. The hnal correction is about, 0.6 for a Conlombic wave function and 

^The squared radial wave function at the origin is of the order of 6 to 6.5 GeV^. 
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Figure 7: Kinematics relevant for the first diagram contributing to the radiative corrections of 
the leptonic decay width. 


0.7 for an Harmonic wave function, relative to the non-relativistic zeroth order calculation 
(Weisskopf-Van Royen approximation). 

These relativistic corrections show that realistic phenomenological two-body potentials, 
such as the Richardson 11 or the Buchmuller-Tye potentials, are able to reproduce the 
decay width of heavy quarkonia, both in the charmonium and the bottomonium sector. This 
solves a long-standing problem in this domain, as phenomenological two-body potentials were 
able to reproduce the whole spectrum of charmonium and bottomonium states, but not their 
leptonic decay width. As already explained in Ref. 0 , the importance of relativistic corrections 
originates from the fact that the leptonic decay widths are very sensitive to the high momentum 
components of the wave function. This, because it involves the integral over all momentum 
range of the wave function and not the square of it, as in many other observables. 
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6 Appendix 

We detail in this appendix the calculation of the radiative corrections to the leptonic decay 
width. Here, we concentrate on the hrst diagram shown in Fig. The second one can be 
calculated analogously. The kinematical conventions are shown in Fig. |^. Applying the dia¬ 
grammatical rules given in Ref. [Q, we hnd: 

= -4, f dr dr' dTp 

V 3 y 2epi(27r)3 2e,.i(27r)3 r r' tq 

^ ^{pI ~ ~ m‘^)5{q^)Q{u!-p2)Q{uj-r2)Q{u!-q){s — M^) 

X + rtPjYifi + m)(j)''{m - f2)'yuim - :^2 + (pr')] , (40) 
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with s = (ri + r 2 )^. The integrations over r, t' and Tq can be done using the delta functions 
induced by the on-mass shell condition of LFD. We also use momentum conservation at vertices, 
including spurion momenta, to express the various momenta. We have: 


with: 


f Xf 2 2'1 

/ — d[r 2 — m ) 

j 

[ ^8{q^) 

J To 


(l-a;)(s-M2) ’ 

1 

{1 - x'){s' - M^) ’ 
1 

2u}-p{x — x')to ’ 


s'={pi +P 2 f , X-- 


OJ-Ti 


s-M^ 

2uj-p 


UJ-p 

s' -M^ 


, X 


U!-pi 


UJ-p 


T 


2uj-p 


, rQ=T 


m — Ti-pi 

uj-p{x — x') 


Thus, we can express s, s' and all other invariant quantities in terms of the integration variables. 
These integration variables can be conveniently chosen as x, x', R R'j, and R t-R't, with: 


R = Ti — xp and R' = pi — x'p , 


(41) 


and Rt is the transverse part of R relatively to the direction of uj, with the consequence 
that R^ = —R^ (see the appendices in Ref. @ for more details). The calculation of the five 
dimensional integral is done numerically using a Monte Carlo procedure. 
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